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Molecular beams of rare gas atoms and D2 have been 
diffracted from 100 nm period SiN x transmission gratings. 
The relative intensities of the diffraction peaks out to the 8th 
order depend on the diffracting particle and are interpreted in 
terms of effective slit widths. These differences have been an- 
alyzed by a new theory which accounts for the long-range van 
der Waals —C3/I 3 interaction of the particles with the walls 
of the grating bars. The values of the C3 constant for two dif- 
ferent gratings are in good agreement and the results exhibit 
the expected linear dependence on the dipole polarizability. 
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Already in 1932 Lennard- Jones [jD predicted that the 
van der Waals interaction of atoms and molecules with 
solid surfaces is given by 



V = -% l>10k 



(1) 



where I is the distance from the surface. This poten- 
tial plays an important role in understanding virtually 
all static (thermodynamical) and dynamical aspects of 
gas adsorption phenomena. Despite its importance, very 
few experimental determinations of C3 have so far been 
reported and most of our present knowledge is based on 
theoretical estimates [0 . The pioneering experiments by 
Raskin and Kusch on the deflection of Cs atoms from a 
conducting metal surface || have recently been extended 
to alkali atoms in high Rydberg states by measuring the 
transmission through 8 mm long narrow (2 — 9 fim) chan- 
nels as a function of their principal quantum number n 
Similar techniques have also been applied to the in- 
teraction of alkali atoms in their ground state J5]]6| or in 
low excited states m . Although the scattering of many 
different atoms and molecules from solid single crystal 
surfaces has been extensively studied, the reflection coef- 
ficients are relatively insensitive to the weak long range 
attractive forces since the collisions are largely deter- 
mined by the reflection from the hard repulsive wall close 
to the surface || . 

Here, a new atom optical technique using transmis- 
sion grating diffraction Pful of molecular beams is em- 
ployed. The van der Waals force causes a change in the 
diffraction intensities just as a smaller slit width would. 
A newly developed theory makes it possible to interpret 



measurements over a range of different beam energies in 
terms of the potential constant C3. For an incident plane 
wave the diffraction peak heights depend on the number 
of illuminated slits N, as iV 2 . With N = 100 slits the 
gain in sensitivity is about four orders of magnitude over 
previous experiments. 

The measurements were made with a previously de- 
scribed 1 10 1 molecular beam diffraction apparatus. The 
beams are produced by a free jet expansion of the puri- 
fied gas through a 5 /im diameter, 2 long orifice from 
a source chamber at a temperature To, into vacuum of 
about 2 x 10~ 4 mbar. At T = 300 K the source pres- 
sure Pq was 140 bar for He, Ne, Ar and D 2 and 50 bar 
for Kr. At lower source temperature P was reduced to 
avoid cluster formation. The atomic beams are charac- 
terized by narrow velocity distributions with Av/v ~ 2.1 
% (He), 5 % (Ne), 7.6 % (D 2 ), 7.7 % (Ar), and 10 % (Kr) 
at To = 300 K, where Av and v denote the full half width 
and the mean value, respectively. After passing through 
the 0.39 mm diameter skimmer the beam is collimated by 
two 10 /jm wide and 5 mm tall slits 6 cm and 48 cm down- 
stream from the source before it impinges on the silicon 
nitride (SiN^) transmission grating with a grating period 
of d = 100 nm and 5 mm high slits with nominal widths 



of s. 



50 nm jLll placed 2.5cm behind the second 



collimating slit. The diffraction pattern is measured by 
rotating the electron impact ionization mass spectrome- 
ter detector around an axis parallel to the grating slits. A 
third, 25 /im wide slit, 52 cm downstream from the grat- 
ing, provides a measured angular resolution of 70 /irad 
(FWHM). 

Transmission measurements with He and Kr atomic 
beams indicate that the grating bars have a truncated 
trapezoidal profile (thickness in the beam direction t) 
|]l2| , ^3| with the narrow face towards the incident beam. 
The measured wedge angles (3 and geometrical slit widths 
so (see below) are listed in Table |. 

The diffraction measurements are illustrated in Fig. |l| 
for four inert gases as a function of the perpendicular 
wave vector transfer k = ksmi9, where -d is the diffrac- 
tion angle. The area under the n-th order diffraction 
peak, J„, is proportional to the grating slit function eval- 
uated at the diffraction angle of the maximum position, 
$ n . For this grating, I, which has equally wide bars and 
slits, the zeros of the slit function coincide with the even 
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diffraction orders H], which are therefore expected to 
vanish. Whereas for He this is almost the case, for the 
heavier rare gases, an increasing deviation is observed. 
For example, the small He intensity ratio of the second 
and third order peaks is slightly larger for Ne, almost 
unity in the case of Ar and, finally, for Kr is greater than 
one. Similar trends are observed for the ratio of the sixth 
and fifth order peaks and in the ratio of the most intense 
zcroth and first orders, which increases significantly from 
about 0.39 for He to about 0.52 for Kr. 

These differences are attributed to the interaction of 
the atoms with the bar walls, Eq. ([!]), which so far has 
not been accounted for in the theory of atom/molecule 
diffraction. For a plane wave e lkz incident on a transmis- 
sion grating with perfectly reflecting grating bars and 
with an additional (attractive) potential at the bar sides, 
the diffracted wave function is, for large r, 



V(r) — /(t?)- 



s i(fcr— 7r/4) 



(2) 



where r 2 = x 2 + z 2 is in the scattering plane normal to 
the height of the slits. The scattering amplitude /(#) is 
determined by the grating transmission function ip(x, 0), 
i. e. by the wave function at the far side slit boundaries 
(z = 0), which depends on the attractive potential. Huy- 
gens' principle E3 yields 



costf 



dx\jj{x 1 0)e 



-ikx sin i5 



(3) 
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If the slit and the bar widths are much larger than the de 
Broglie wave length A, the intensity 1(d) — \f('d)\ 2 can 
be written as a product 

2 



(\Nkdsm.d) 



sin (5 fee? sin $J 



i/ s iitwr 



(4) 



where N denotes the number of slits and |/ s iit| 2 is the slit 
function. Thus, the atomic diffraction pattern consists 
of principal maxima at the diffraction angles sin?9„ = 
nX/d, n — 0, ±1,±2, ... while |/ s nt(i?)| 2 plays the role 
of an envelope function. Eq. gives, after a change of 
variable from x to a variable with the origin at the edge 
of a slit, Q = sq/2 — x, 



ji9 



d( cos 



/so 



c 



r(C), (5) 



where r(C) = ^(V 2 - C, 0), < ( < s /2, is the single- 
slit transmission function. 

It is instructive to first deduce the general structural 
form of /siit (I?)- Since the grating bars reflect those atoms 
which touch the bar walls, the wave function in the slit 
vanishes at the walls, i. e. r(0) = 0. Taking this into 
account and after a partial integration Eq. (^J) becomes 



/slitW = 



cosfi 



(I) 



(6) 



where 

with $(0) = 1. The logarithm of <E> can be expanded as 



(7) 



OO r . . \ 



(8) 



where the complex R n are known as cumulants |lf 



Ri = 



d((- 



r/ (C) 
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(9) 



etc.. For the small wave- vector transfer n of interest here, 
only the first two terms are needed in the series Eq. (||). 
The single-slit amplitude Eq. (H) then becomes 



„cosi? 

/ S lit0?) = 



(?) 



ill)] 



(10) 



For a comparison with experiment the surface rough- 
ness of the grating bars must be accounted for. In a 
first approximation roughness has been included by rigid 
shifts of the individual bar sides (see also Ref. p6|), 
which are randomly Gaussian distributed. In the case 
of a weak surface potential, this results in an additional 
Debye- Waller like damping factor exp(— fc 2 <r 2 sin 2 i? n ) in 
the intensity ratio of the principal maxima, I n / Iq, where 
Co is the variance of the geometrical slit width (T^| . Tak- 
ing this into account, Eq. (0) with Eq. (fof) yields 
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where a = CTq + Re(i? 2 ), s c ff = s o — 2Re(i?j) and 5 = 
2Im(i?i). The first term in the brackets of Eq. (nW leads 
to a Kirchhoff-like slit function (see e. g. Ref. [^3|) with 
a Debye- Waller term and an effective reduced slit width 
s e ff, while the second term suppresses the zeros of the 
Kirchhoff pattern, as can be seen in the insets of Fig. |]. 

The effective variance a 2 as well as s e g and i5 in 
Eq. ( |ll| ) can be calculated for the potential Eq. (Q). 
The standard eikonal approximation 17 l| is used to 
determine the grating transmission function, given by 



V>(x,0) 



,iip{x) m g]|^ s anc j zer0 e i sew h ere _ The 



phase shift reads 



(p(x) 



1 

hv 



dzV(x, z), 



(12) 



where v = Hk/m is the particle velocity. Taking the 
trapezoidal bar profile into account, after some algebra 
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the single-slit transmission function becomes 



r(C) = exp 



t COS P C3 



hv 



1 - 



5^ tan/3 

T2 

I tan (3 J 



(13) 



An analysis of Eqs. ( [13| ) and (g) reveals that Re(i?i) and 
hence s e ff is especially sensitive to the potential. 

The effective slit width s e ff as well as 8 and a were 
determined from the experiment by fitting the relative 
experimental diffraction intensities I n / 1\ as depicted in 
the insets of Fig. [I] to the corresponding ratios determined 
from Eq. (|ll|). These ratios and not I n /Io are compared 
with theory since small concentrations of clusters in the 
beams can falsify the Iq intensities. The effective slit 
widths are plotted versus the particle velocity in Fig. || 
(points) for two different gratings. The difference be- 
tween the effective slit widths for To = 300 K beams and 
the geometrical slit width sq increases from 1 nm (He) 
to more than 6 nm for Kr as expected from the increas- 
ing interaction strength of the van der Waals potential. 
With increasing C3 the slope of the curves also increases. 
The solid lines in Fig. || represent least squares fits of the 
theoretical expression s e s — sq — 2Re(i?i), with Ri given 
by Eqs. @ and ([b|), to the experimentally determined 
effective slit widths, which allow for the determination of 
C3 and So- Since He has the smallest polarizability and 
measurements over the largest range of velocities were 
possible they were used to determine the values of sq in 
Table | for each of the gratings. Identical values for sq 
were obtained from D2 measurements. This value of sq 
was then used for Ne, Ar and Kr, with C3 the only re- 
maining fit parameter, and hence for these systems mea- 
surements at various velocities are not necessary. 

The C3 parameters are plotted versus the static elec- 
tric dipole polarizabilities a in Fig. [| The error bars 
were determined by assuming a realistic uncertainty in 
the bar geometry by varying (3 by ±2° in Eq. (13). This 
uncertainty seems to be the only systematic source of er- 
ror in the present C3 determination and leads to errors of 
about 20 %. Figure || indicates that the influence of the 
surface potential is restricted to distances much smaller 
than the slit width and therefore, by Ref. p8[ , corrections 
due to the finite bar width should be negligible. 

Within the error bars the data from both gratings fall 
on a straight line in agreement with Hoinkes' empirical 
rule jq]. Accordingly the slope provides information on 
the optical dielectric constant of the grating material. 
An approximation to the theoretical expression for C3 
predicts that D2 should in fact have a slightly smaller 
ratio of C%/a than the rare gas atoms, while among them 
Ne is expected to have the largest ratio. It is satisfying 
to see that the small deviations from the straight line in 
Fig. agree with this expected trend. 

The big advantage of the present method is its large 
sensitivity as can be seen from Fig. and its universal- 



ity. In principle all atoms and molecules are accessible 
for study. The only restrictions will be to produce grat- 
ings of different solids and molecular beams with suffi- 
ciently narrow velocity distributions and to reduce the 
corresponding background in the mass spectrometer de- 
tector to assure an adequate signal to noise ratio. The 
present work also allows for a quantitative understanding 
of diffraction intensities in atom optics and atom inter- 
ferometry experiments using transmission structures as 
optical elements. 

We are greatly indebted to Tim Savas and Henry I. 
Smith of MIT for providing the transmission gratings to 
us. Further, we thank Dick Manson and G. Schmahl for 
fruitful discussions. 
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FIG. 1. Diffraction patterns measured with Grating I for 
He, Ne, Ar, and Kr at the same beam energy (To = 300 K). 
The insets contain a comparison between least-squares fits of 
I n /Ii determined from Eq. ( p"l| ) with continuous values of n 
(solid lines) and Kirchhoff theory (dashed lines) to measured 
diffraction intensity ratios (points). 



FIG. 2. Effective slit widths plotted as a function of the 
particle velocity for He, Ne, D2, Ar, and Kr beams. The solid 
lines are theoretical curves determined from Eqs. ([]) and ( |l3| ) 
with the C3 parameters in Figj3[ Data points indicate fits 
of I n /h determined from Eq. ( |1 1| ) to experimental intensity 
ratios obtained from diffraction measurements with two grat- 
ings. 

FIG. 3. Measured C3 values of silicon nitride (SiN^) ob- 
tained in this work plotted as a function of the static electric 
dipole polarizability of the respective atom, a (see Ref. [^]). 
The solid line is a linear fit of the data. 



TABLE I. Geometrical properties of the three gratings. 



grating 


P [°] 


so [nm] 


I 

II 

III 


7.5T2 
8.7T2 
12.7T2 


50 
67.5T0.1 
71.2±0.1 



4 




Wave Vector Transfer k = /csind [ 2tz / d ] 



